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Abstract: Magnetism is a ubiquitous property of astrophysical plasmas, yet stellar magnetism still remains 
far from being completely understood. In this review, we describe recent observational and modelling 
efforts and progress to expand our knowledge of the magnetic properties of high-mass stars. Several 
mechanisms (magneto-convection, mass-loss quenching, internal angular momentum transport, and 
magnetic braking) have significant implications for stellar evolution, populations, and end-products. 
Consequently, it remains an urgent issue to address and resolve open questions related to magnetism in 
high-mass stars. 


Keywords: magnetism; high-mass stars; stellar evolution; magnetic field evolution; fossil fields; dynamo 
fields; angular momentum transport; magnetic braking; mass-loss quenching 


1. Introduction 


Magnetism plays an important role in several astrophysical phenomena. For example, 
stellar and planetary magnetic fields may be a crucial aspect of planetary habitability, e.g., [1,2]. 
Magnetic fields may significantly impact star and planet formation, e.g., [3,4]. The magnetic 
characteristics of neutron stars have been extensively studied via observations, e.g., [5,6] 
and numerical approaches, e.g., [7,8], and the magnetic properties of the interstellar medium, 
galaxies, and galaxy clusters have also been subject to ongoing investigations, e.g., [9]. However, 
the magnetic properties of non-degenerate stars are still not fully understood. Cowling and 
Spitzer concluded that, based on the available magnetic flux during star formation, it should 
be possible that some remnant, so-called fossil fields, remains in stars [10,11]. Such fossil 
fields will be discussed in detail in this review. Moreover, the dynamical characteristics of 
stars also make them plausible sites to generate in situ magnetic fields. A contemporaneously 
generated (and re-generated) magnetic field taps from a kinetic energy reservoir and converts it 
to magnetic energy. This is the fundamental principle of a dynamo. Various dynamo processes 
have also been proposed to operate in stars, and we will address this topic too. In this review, 
we will focus on high-mass stars with O and B spectral types. They are characterised by a 
convective core and a radiative envelope on the main sequence. This structural property is 
fundamental to unveil their magnetic characteristics. High-mass stars play a vital role in the 
chemical evolution of the Universe, e.g., [12]. Via their powerful stellar winds and ionising 
radiation, they regulate star and planet formation by stellar feedback, e.g., [13,14]. Stars initially 
more massive than about 8 Mo leave behind neutron stars and black holes by the end of their 
evolution. Understanding the structure and evolution of the progenitors of compact remnants 
is required to explain several astrophysical phenomena, which range from observations of 
fast-radio bursts, e.g., [15-17] to detections of gravitational wave events, e.g., [18]. 
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2. Observations 


In this section, we briefly address some milestones from observational studies as well 
as important new clues that help us shape theoretical advancements. For a comprehensive 
background on observational aspects applicable to various kinds of stars, including the Zeeman 
effect, polarimetry, and modern spectral convolution and imaging techniques, we refer the 
reader to [19-21]. 


2.1. Historical Perspective 


The first confirmation that magnetism is present in stars took place shortly after Zeeman’s 
discovery of spectral line splitting under external magnetic fields. Hale could measure the line 
splitting in sunspots and infer their magnetic field strengths [22]. Since then, our knowledge of 
solar activity has significantly improved, and it was understood that the underlying origin of 
the observed magnetism is related to convection and rotation, e.g., [23,24], although the exact 
role of convection is still debated, e.g., [25]. In fact, while there still remain open questions 
about the solar cycle and the solar magnetic field generation and activity, given the supporting 
observational evidence, it is now generally accepted that all solar-type stars (and, more gen- 
erally, stars with convective surfaces, including M dwarfs, Red Giants, and Red Supergiants) 
produce a dynamo mechanism given the ingredients of convection and differential rotation, 
e.g., [21,26,27]. This significant finding further raises the question of whether any convectively 
unstable region in stars could operate a dynamo cycle and thus produce magnetic activity 
hidden from the observer. 


2.2. From Ap/Bp to OBA Stars 


The first magnetic field on a star other than the Sun was detected by Babcock on a chemi- 
cally peculiar star, 78 Vir [28]. In subsequent years, the discovery of chemically-peculiar Ap/Bp 
stars grew further, e.g., [29,30]. It remains a very rich field of study, with increasingly more 
sophisticated techniques to infer surface spots and abundance patterns, e.g., [31-35]. There is 
a demonstrable link between the magnetic characteristics of magnetised pre-main sequence 
Herbig Ae/Be stars and main sequence Ap/Bp stars [36-39], suggesting that the observed 
magnetic fields of main sequence stars are already developed in the pre-main sequence phase. 
This observational evidence is a crucial piece of information in ongoing investigations about 
the origin of fossil fields. 

One indirect hint that other early-type main-sequence stars with radiative envelopes 
may harbour surface magnetic fields originated from the X-ray properties and evidence of 
material falling back onto the stellar surface of ©! Ori C, e.g., [40-42]. Analytical wind-shock 
models were developed in the context of a stellar magnetosphere [43]. These predictions were 
also verified by state-of-the-art simulations studying the interaction of stellar winds and the 
magnetosphere [44], as well as further X-ray observations [45]. The X-ray characterisation of 
magnetic massive stars remains invaluable, e.g., [46—48]. Importantly, these predictions also 
led to optical spectropolarimetric observations (using circular polarisation, that is, the Stokes 
V parameter), which detected a surface magnetic field in ©! Ori C, a chemically non-peculiar 
hot star [49]. Since then, large-scale spectropolarimetric surveys have been conducted to study 
the magnetic properties of stars with O, B, and A spectral types [50-59]. This led to valuable 
statistical information on magnetism in high-mass stars. The surveys concluded that magnetism 
is unambiguously detected in about 10 percent of the sample stars, regardless of spectral type. 
The main common ingredient seems to be the similarity in stellar structure (convective core 
and a radiative envelope), but dependence on stellar parameters (mass, rotation) could not 
be found [59,60]. It has also become clear that chemical peculiarities may not be evidenced 
in some magnetic stars. The lack of such peculiarities is generally attributed to stellar winds, 
which can prevent vertical chemical stratification in the upper stellar layers. Indeed, many 
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magnetic O and B-type stars do not exhibit the same kind of peculiarities as Ap and Bp stars, 
yet their magnetic characteristics are broadly consistent. 

The empirical characterisation of these magnetic fields suggests that the typical structure is 
predominantly an oblique dipole (misaligned magnetic and rotation axes) in the majority of the 
cases, although extensive spectroscopic and spectropolarimetric monitoring may help reveal 
departures from the first-order pure dipole fits, e.g., [61]. An important clue about the possible 
origin of magnetism in high-mass stars relates to recent measurements of obliquity angles 
by [62]. The findings indicate a random distribution in the tilt of the magnetic axis compared 
to the rotation axis of the star. This also indicates that magnetic precession may have a longer 
timescale than the main sequence evolution of the star [63,64]. In addition, a small fraction 
of known magnetic hot stars cannot be reconciled with dominantly dipolar geometries: their 
peak magnetic field strengths are in higher-order spherical harmonic components, suggesting 
a complex magnetic geometry, e.g., [65,66]. Interestingly, these stars are generally consistent 
with a young stellar age. 

The inferred magnetic field strengths range from a few hundred G to up to tens of kG. The 
strongest magnetic field reported (20 kG inferred dipolar strength) for a seemingly chemically 
normal O-type star is of NGC 1624-2 [61,67]. The detection of “ultra-weak” magnetic fields in 
Vega [68,69], an intermediate-mass A-type star, has been surprising. A few other A-type stars 
have been discovered to host ultra-weak fields, although they possess chemical peculiarities 
typically associated with the Am class [70,71]. It is unknown if such weak magnetism may or 
may not be a widespread phenomenon since the origin of such fields is unclear. However, the 
observed strong fields are now commonly explained as fossil remnants [72-75]. The bi-modality 
of observed magnetic field strengths is separated by a “magnetic desert”, i.e., a scarcity of 
inferred dipolar magnetic field strengths in the range from a few G to about three hundred 
G [76,77]. It seems that this cutoff separates magnetic fields that may either be produced by 
different mechanisms or may reflect on the distribution of initial magnetic field strengths [78]. 

The results from spectropolarimetric surveys raise important questions that are major 
drivers for current investigations: Why do only a fraction of hot stars show evidence of strong 
magnetic fields? How do these fields evolve over time? How does the magnetic geometry 
change over time? 

To gain further insights from observations, it is essential to increase the sample size of 
known magnetic high-mass stars; therefore, high-yield surveys are much desired. A successful 
strategy has been the identification of characteristic light curve variations in photometric data 
and subsequent spectropolarimetric follow-up observations, e.g., [79-82]. Recent advance- 
ments include further diagnostics from infrared observations , e.g., [83-85], radio observations, 
e.g., [86-90], and planned UV spectropolarimetric missions, e.g., [91-93]. The multiwavelength 
characterisation is invaluable to fully explore the physical behaviour of observed magnetic 
fields. One incredibly important piece of missing information, the metallicity dependence of 
observed magnetism [94,95], is beyond reach with current instrumentation. High-resolution op- 
tical spectropolarimetric units, suitable for stellar magnetometry of hot stars, are only mounted 
on 4 m class telescopes. To detect magnetic high-mass stars in the Magellanic Clouds and 
beyond, it is essential to include high-resolution spectropolarimetric units on 10 m class and 
larger aperture telescopes. 


2.3. Neutron Stars and Magnetars 


Although the detailed discussion of compact objects is beyond the scope of this review, 
we should outline that there is ongoing debate as to how the magnetic field of these objects 
is acquired. In particular, a subclass of neutron stars, magnetars, shows evidence of surface 
magnetic fields of the order of 101-10!" G, whereas the magnetic field strengths of typical 
neutron stars range between 10!! and 10! G [5]. 
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One scenario suggests that the magnetic flux of main sequence stars is conserved and 
serves directly as the magnetic field of the remaining compact object [75,96,97]. A major 
challenge for this scenario is that magnetic flux should be essentially entirely preserved during 
stellar evolution and core collapse. It is unclear if these assumptions would hold. Another 
issue is that the fraction of strongly magnetised main sequence stars may not match the fraction 
of strongly magnetised neutron stars [98]. 

Other origin scenarios involve in situ amplification from a seed magnetic field at around 
the time of core collapse via a dynamo mechanism [99]. This could be realised via a convective 
dynamo, e.g., [100], a magneto-rotational instability-driven mechanism, e.g., [101-103], or even 
a combination of these [104]. New approaches have also considered a Tayler-Spruit mechanism 
to further amplify a 101? G magnetic field in a few seconds to 1016 G [105]. 

In order to understand the origin of strong magnetic fields in magnetars, it remains 
an essential task to map out the magnetic characteristics of their progenitors and constrain 
how much magnetic flux may be available at the time of core collapse. This could also have 
implications for detections of fast-radio bursts that have been linked to magnetars, e.g., [15-17]. 


3. Theoretical Background 


In this section, we overview some key concepts that form the foundation of new modelling 
approaches. These concern the kinds, stability, and effects of magnetic fields. There are several 
excellent textbooks and reviews on the topic—for example, [106-112]. 


3.1. Fossil vs. Dynamo Fields 


Stars could inherit some magnetic flux from the molecular cloud. The parent molecular 
cloud has its own magnetic field. As the cloud collapses, the magnetic field lines are compressed, 
and the field strength increases. This amplification is a consequence of the “frozen-in” flux 
condition (see below). In fact, it could be shown that some magnetic flux must be dissipated (via 
non-ideal MHD processes) during this phase to avoid a “flux catastrophe” problem of ending 
up with too much magnetic flux in stars [3,113]. A remnant magnetic field from star formation 
can relax into an equilibrium state and form a large-scale structure. Demonstrating the stability 
of such fields has been a pioneering work by Braithwaite and collaborators [73,114-116] and 
has also been recently revisited [117,118]. 

Magnetic fields induced by a dynamo cycle are fundamentally different from the above 
picture. A dynamo cycle converts kinetic energy to magnetic energy and thus operates on a 
completely different time scale. Dynamo theory has been extensively developed; for a com- 
prehensive review, see, e.g., [108]. Several dynamo mechanisms describe a local phenomenon. 
Therefore, dynamo-produced fields are often characterised by length scales that are relevant to 
the powering mechanism (e.g., convection). The perturbed magnetic fields can be considered as 
small-scale, for example, in comparison to the size of the star. An electromotive force is necessary 
to reconstruct a large-scale field from small-scale perturbations. This is often described via the 
mean-field theory [106]. There are nonetheless still debates about relevant length scales and as 
such whether some dynamos would be small- or large-scale. 


3.2. Time Evolution of Magnetic Fields 


A general MHD equation of the time evolution of magnetic fields can be cast in the form 


of: 
ƏB 
57 = V x(u xB+1V xB), (1) 


where B is the magnetic field vector, t is the time, V is the spatial derivative in three dimensions, 
u is the fluid velocity, and y is the magnetic diffusivity. 
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3.2.1. Fossil Fields 


In the case of relaxed, equilibrium, fossil magnetic fields, the above equation simplifies to: 


& =nV°B, (2) 
assuming that the magnetic diffusivity has no spatial dependence. There is no active generation 
or induction term. Assuming that the most relevant physical process leading to the dissipation 
of fossil fields is Ohmic diffusion, the evolution of fossil fields is expected on a very slow 
(Ohmic) time scale. It can be shown that, for the Sun, a general Ohmic dissipation time scale 
(Tonm © 17/7, with I being a characteristics length scale) is longer than the nuclear-burning 
time scale [10]. The accurate calculation of Ohmic dissipation is difficult in the context of 
high-mass stars. The magnetic diffusivity varies as a function of radial distance inside the star 
and depends on thermodynamic quantities. For example, the temperature profile drops several 
orders of magnitude from the core to the surface. Thus, current estimates for the magnetic 
diffusivity remain quite uncertain [119]. It is necessary to understand how the fossil field is 
anchored inside the star to obtain better quantitative estimates of dissipation processes. In 
high-mass stars, the nuclear-burning time scale is much shorter than in the Sun (~ few Myr vs. 
~10 Gyr). If the Ohmic dissipation time scale for the Sun and high-mass stars were roughly 
similar, any fossil field in high-mass stars should be long-lasting. Some additional mechanisms 
or phenomena still cannot be ruled out as potential sources of magnetic field dissipation. An 
interesting possibility may be tidal dissipation by a companion star [120]. However, rapid 
dissipation requires a close companion on a short-period (<20 d) orbit. Although multiplicity 
is common amongst high-mass stars, the number of short-period systems is a more moderate 
fraction. Thus, it is not immediately evident what fraction of stars are affected by this process. 
In general, the dissipation of fossil fields is expected to be slower than the star’s evolution. The 
principle of the frozen-in flux condition (Alfvén’s theorem [121]) can be applied to deduce the 
surface magnetic field strength. The magnetic flux over a closed surface area S is defined as: 


p= f BdS ~ 47R2B, (3) 
S 


where the unsigned magnetic flux is considered for the surface of a sphere with radius R. 
Since the time derivative of the magnetic flux is zero, the above expression must be constant, 
requiring a relation between the magnetic field strength and stellar radius. If we define the 
surface strength of the magnetic field in a time-dependent manner, we can write: 


Reare (to) 

Bsurf(t) = Beurt(to) 22 — 4 

surf ( ) surf ( 0) R? a) (4) 
with Rsurf being the stellar radius and tọ an initial time (typically assumed to be the Zero Age 
Main Sequence). Since the radius of the star changes during its evolution, we may anticipate a 
change in the surface field strength in response. Some observational studies are consistent with 
this picture [56,58,122], while others infer a more rapid decline of the magnetic field strength 
by a magnetic flux decay [54,59]. These approaches often rely on stellar ages inferred from 
non-magnetic evolutionary models. This source of uncertainty can now be improved upon 
with models incorporating fossil field effects [123,124]. 
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3.2.2. Dynamo Fields 


Dynamo fields have been proposed in several forms [107]. For example, the induction 
equation for a convection and rotation-driven, so-called a-CQ dynamo can be written as: 


— =Vx(uxB+aB), (5) 


where, for simplicity, we assumed no dissipative terms. Such dynamos are mostly constrained 
by the Rossby number Ro = 27t/OTconyv, with angular velocity Q and the convective turnover 
time scale defined as the ratio of a characteristic length scale and convective fluid velocity 
Teonv = l/Veonv, €.g., recently [125]. This implies that such dynamos will show a dependence on 
the rotation rate of the star. In Equation (5), the « effect is responsible for turbulent convection, 
which displaces a toroidal magnetic field and converts it into a poloidal structure. The O 
term (encompassed in the fluid velocity vector u) accounts for shearing the magnetic field by 
differential rotation, which regenerates a toroidal component from a poloidal one. Variations of 
this kind of dynamo include so-called «2 — Q and g? dynamos [107]. The evolution of dynamo 
fields takes place on a short, Alfvén time scale, which is characterised by tg ~ 1/v,4 with 
vA = Brms/\/479 the Alfvén speed. Typical time scales in stars are of the order of a few days 
to a few decades. 


3.3. Stability Criteria for Fossil Fields 


While equilibrium states were traditionally found in several forms, the stability of fossil 
fields has been more challenging. Rigorous studies of the stability of fossil fields have been 
conducted in numerical simulations [114-118]. Several stable configurations exist. Figure 1 
shows two of them: an axisymmetric dipole and a non-axisymmetric, more complex magnetic 
structure [115,117]. Some stable configurations also exist without any observable surface 
manifestation. 


Figure 1. Examples of stable fossil field configurations. Left: Axisymmetric dipole field [117]. The fossil 
field is in stable equilibrium after more than 600 Alfvén times. Coloured lines represent magnetic field 
lines. The grey colour corresponds to the stellar photosphere, and the magnetic axis is shown with the 
red arrow. Right: Non-axisymmetric, complex magnetic structure in a stable configuration after several 
tens of Alfvén times [115]. The blue background represents the star as a spherical gas. Red lines indicate 
magnetic field lines. 
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Figure 2. Magnetic energy as a fraction of its initial value over Alfvén time scales. The polytropic index 
y = dln P/dIno is varied. In a barotropic case, it equals the adiabatic exponent, y = IT = 5/3. The 
magnetic energy declines without reaching a stable configuration in this case, from [118]. 


These studies confirmed early analytical constraints that both purely toroidal [126,127] and 
purely poloidal [128-130] magnetic field configurations are unstable. Following the analytical 
model of Prendergast [131], a stable configuration in numerical models was found by requiring 
a twisted torus shape of mixed toroidal and poloidal modes, in which the toroidal field is 
confined by poloidal lines [73]. In barotropic cases, the magnetic field can be calculated from a 
Grad-Shafranov type equation [132,133]. This concerns the stellar interiors and has also been 
found to be an equilibrium configuration in the presence of rotation [134,135]. However, the 
stability of Prendergast-type fields over the course of stellar evolution has been questioned 
recently [136]. 

Some controversy also concerns the use of a barotropic equation of state. Using analytical 
works and numerical simulations, it has been argued that non-barotropic conditions may be 
important for stable stratification and thus for the stability of the magnetic field; see Figure 2 [7, 
118,137-139]. Barotropic refers to regions where the density depends only on the pressure 
(oe = p(P), or in other words, isobars coincide with iso-density regions) but not on other 
thermodynamic quantities. 

For a realistic radiative stellar envelope, stable stratification, an outward decreasing 
density and an increasing specific entropy profile are necessary. According to [7], this does 
require non-barotropic conditions, i.e., knowledge of other thermodynamic quantities, for 
example, the temperature or entropy in order to remain stable against convection. 

The dipole mode may form a magnetosphere around the star; however, any toroidal field 
diffuses above the stellar surface such that the field can eventually become force free [111,140], 
meaning that the Lorentz force is zero, (V x B) x B = 0. While the magnetic field above the 
stellar surface can possess an azimuthal component (for example, due to stellar rotation shearing 
a pure dipole field), enclosed toroidal fields above the stellar surface are unstable and rapidly 
diffusive into the poloidal component. 


3.3.1. Momentum Equation 
Let us write the force balance, the MHD Euler equation of motion, as: 


v= VP +pg+z-(V x B)xB+20xu (6) 
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where P is the pressure, g is the gravitational acceleration, B is the magnetic field vector, Q is 
the angular velocity, and u is the fluid velocity. Inside stars, the magnetic field is non-force free, 
and there is a non-zero Lorentz force. 

A fossil magnetic field must be stable against various instabilities; otherwise, it would 
decay and become unobservable. Here, we describe two interactions that are key to the stability 
of fossil magnetic fields. In both cases, the stability of the fossil field relies on a critical value. 
It is thus required for an initial magnetic field to have a critical strength in order to “survive” 
these interactions and relax into a stable equilibrium configuration. 


3.3.2. Convective Expulsion and Inhibition of Convection 


A stability criterion for magneto-convection, the interplay between convection and mag- 
netism, was first given by [141]. A first-order estimate can be posited by considering that 
convection in stellar layers can be suppressed by a critical magnetic field when the magnetic 
pressure reaches the thermal pressure: 


ab? = g>? (7) 


that is, when the magnetic field has equipartition strength, e.g., [76,142]. Recently, a more 
sophisticated model was considered, including radiation pressure [143]. The revised stability 
condition is formulated as: 


2 
vA i dinT, 


where Q is the thermal expansion coefficient, V and Vq are the structural and adiabatic 
gradients, v4 and cs are the Alfvén and sound speeds, T4 is the first adiabatic exponent, and 
P is the local pressure. A minimum Alfvén speed v4 = B/,/47¢p, or equivalently a critical 
magnetic field strength, is needed for the above relation to be satisfied. If B > Berit, convection 
can be suppressed or even inhibited, e.g., [144]. This is commonly considered in the context 
of dynamo fields in convective regions of low-mass stars (e.g., [145]). On the other hand, if 
B < Bait, convection will expel the magnetic field outside of the convectively unstable region. 
This is commonly known as convective expulsion, e.g., [146-148]. 


3.3.3. Rotational Expulsion and Inhibition of Differential Rotation 


Similar to the above constraints, differential rotation could also lead to the local removal 
of an initial magnetic field if the difference in rotation frequency between adjacent layers 
is much larger than the Alfvén frequency [111,149]. In this case, the magnetic field may be 
wound-up such that a high magnetic diffusion effectively dissipates it, leaving no magnetic 
flux in differentially-rotating regions. 

Once a fossil field is formed, the rotation rate of the star will slow due to magnetic 
braking. Importantly, strong differential rotation cannot develop in those regions where the 
field spreads [106,150]. Instead, the radial rotation profile becomes uniform [151]. Therefore, 
models with fossil fields are unlikely to be compatible with strong differential rotation, unless 
they exist in different regions of a star. 


3.4. Angular Momentum Transport and Loss 


Magnetic field lines are very efficient at transporting angular momentum. We can define 

a magnetic stress tensor (Maxwell stress tensor) T as a quantity whose negative divergence 
yields the Lorentz force: 

FLorentz =-V.-T. (9) 
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The Maxwell stress tensor leads to a momentum flux (force per unit area), also known as 
Poynting flux. The r, @ component of the stress tensor is given by: 
ByBy 
Tro = j 
"? ~~ “Ar 


(10) 


with B, and Bọ the radial and azimuthal component of the magnetic field vector, respectively. 
Integrating this expression over a solid angle yields the azimuthal magnetic force. The scalar 
value of the magnetic torque is then: 


N = B,Bgr?rsin (11) 


with r the radius and 0 the co-latitude. This is equivalent to the rate of change in angular 
momentum per unit time. 

There are different scaling relations in the literature to express B, and By, depending on 
the problem (for example, the saturation amplitude of a dynamo field) or simply the kind of 
magnetic field considered. Consequently, expressions for the magnetic viscosity (or resulting 
magnetic diffusion) can vary, e.g., [25,152-154]. In turn, an exact magnetic angular momentum 
transport equation, which is usually assumed to be a diffusive process, is still somewhat 
uncertain [155]. It has been argued that the magnetic Maxwell stress leads to an advective term 
in the angular momentum transport equation, which could aim to drive the system towards 
uniform specific angular momentum instead of uniform angular velocity [156]. The diffusivity, 
whether produced directly via Maxwell stresses or indirectly via magnetic turbulence, must 
be high. Indeed, it is expected that poloidal field lines could “freeze” rotation and enforce 
solid-body rotation (Section 3.3.3). In addition, magnetic fields are also expected to affect wave 
propagation in stellar interiors, e.g., [157-159]. 

Considering that there is an outward mass flux from stars, Weber and Davis [160] showed 
that the net rate of change in angular momentum per unit time caused by the solar wind can 
be cast in the form of: a2 

; 2 

an 3MO.R i (12) 
where M is the mass-loss rate, O, is the surface angular velocity, R4 is the Alfvén radius, and 
the factor 2/3 results from integrating over the sphere. The occurring mass loss implies angular 
momentum loss, thus the above formula is commonly adopted to describe the stellar spin down 
via magnetic braking [161,162]. In the Weber and Davis model, R4 is defined as a location 
where the Alfvénic Mach number is unity. In pioneering MHD simulations of massive stars, 
it was shown that this analytical scaling relation remains a good approximation for the net 
rate of angular momentum loss, albeit changing the expression for the Alfvén radius from the 
split monopole model to a more appropriate dipole (or more general, higher-order harmonic) 
scaling [163,164]. Subsequent improvements concern including the tilt angle of the magnetic 
field, which only modestly impacts the above scaling relation [165]. 


4. Modelling Approaches 
4.1. Dynamos in the Stellar Core 


For the example of solar-type magnetism, it is nowadays commonly expected that the 
convective cores of massive stars are good candidates to induce a similar « — OQ dynamo cycle 
and produce magnetic fields. Such convection and rotation-driven dynamos were first studied 
via semi-analytical methods and numerical one-dimensional stellar evolution models [119]. 
A central question has been whether magnetic fields produced in the core could reach the 
stellar surface. However, opposite conclusions were reached regarding the buoyant rise time of 
magnetic flux tubes, and this question remains unsettled [166,167]. 
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Thus far, a few numerical MHD simulations have targeted this domain. The convective 
cores of A and B-type star models were studied, revealing strong dynamo activity [168-170]. 
The first MHD model of a near-core magnetic field of a B-type star, supported by asteroseismic 
inference, has recently been studied by [159]. The inferred field strength is unlikely to be 
reached by a fossil field, supporting a convective dynamo origin. 


4.2. Dynamos in Convective Shells 


In recent years, it has become possible to perform multidimensional MHD simulations 
focusing on silicon/oxygen shell burning in massive stars [171]. This stage is within a few 
minutes of the star reaching core collapse. The rotational and magnetic properties at this 
stage can have a large impact on the core dynamics. For example, asphericities may aid the 
neutrino-driven explosion. The magnetic field might also become the seed of a subsequent 
dynamo action to create the neutron star’s magnetic field. The inclusion of magnetic fields is a 
great advancement in such models in addition to previous hydrodynamic simulations [172- 
175]. Such MHD simulations are computationally expensive and still require a comprehensive 
approach to consider several physical ingredients, for example, stellar rotation. Nonetheless, 
their impact on compact object formation may have far-reaching consequences. 


4.3. Dynamos in the Stably Stratified Envelope 


The envelopes of high-mass stars on the main sequence are dominated by radiative energy 
transport. In these stably stratified regions, turbulent convection is absent, hence it cannot 
support the a mechanism to operate a dynamo. However, there may exist other ways of 
supporting a dynamo cycle in radiative stellar regions. A breakthrough in this field took 
place when Spruit proposed a dynamo that operates via the Tayler instability and differential 
rotation [25,126]. Spruit argues that in general the a effect of convection may be replaced by 
instabilities in the magnetic field itself [25]. This requires that differential rotation, quantified in 
the form of q = dInQ/dInr, must reach a sufficient value. In this picture, it is considered that 
the first instability to set in is the Tayler instability, and it effectively takes the role of turbulent 
convection as the « mechanism to produce a contemporaneous dynamo cycle [25]. 


4.3.1. Tayler Instability 


The basis of the Tayler instability, also known as Pitts—Tayler instability, is a pinch type 
deformation of the magnetic field [126,130,176]. It was found that a purely toroidal field 
becomes unstable when reaching a threshold value and starts to produce perturbed fields in the 
radial and latitudinal directions, that is, diffusing into small-scale poloidal field components. 
The Tayler instability is a fundamental piece to understand a possible dynamo cycle, yet it 
has only recently received some more attention via using analytical methods [177] and 3D 
simulations [178]. 

Spruit proposed that a dynamo mechanism could operate by i) differential rotation wind- 
ing up a poloidal field and generating a toroidal field, and ii) the Tayler instability re-generating 
the poloidal field [25]. It is necessary that the small-scale perturbations are amplified via an 
electromotive force (u x B) to re-create the poloidal field. The key point to obtain a stable 
dynamo cycle is that the amplification must be balanced by some form of dissipation. Amongst 
other things, this remains a subject of debate. In the work of [25], the dissipation concerns a 
thermal diffusion of the large-scale toroidal field to saturate the dynamo, and numerical models 
of [114] claimed to confirm this cycle. 

However, subsequent works have questioned whether the dynamo cycle could operate. In 
other numerical simulations, it was found that the dynamo cycle does not reach a steady-state 
since the electromotive force is insufficient to regenerate the large-scale poloidal field [179]. It 
was also pointed out that the Tayler instability generates local perturbations; thus, a large-scale 
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dynamo as envisioned, on the length scale of the stellar radius, may not be achieved [180]. The 
issue of the length scale on which the Tayler instability operates is still not settled [155]. 
Recently, adjustments were made to Spruit’s model by supposing that the dissipation 
concerns the small-scale perturbed fields and not primarily the large-scale toroidal field [153]. 
This way, higher saturation strengths can be reached and consequently the resulting magnetic 
viscosity is also higher. This allows for spinning down the stellar cores to quite low rates. 


4.3.2. Magneto-Rotational Instability 


Another branch of possible dynamo solutions is related to the magneto-rotational instability 
(MRI), which was introduced by [181] for accretion disks following earlier work of [182]. It has 
also been studied in core-collapse supernova models, e.g., [101,103,104,183,184]. Nevertheless, 
the consideration of this mechanism in main-sequence stellar models has only recently been 
gaining more momentum [185-188]. The key principle is that differential rotation must reach 
a sufficient strength to trigger the instability. A negative angular velocity gradient is required 
such that the angular velocity increases towards the centre of the star. The MRI is a very rapid 
exponential magnetic field amplification process that primarily generates a strong toroidal field. 
It remains a matter of debate whether in realistic stellar conditions the Tayler instability or the 
magneto-rotational instability could occur first. 


4.3.3. Mean-Field Dynamos 


Further advances in this field concern dynamo models that differ more substantially from 
MRI or Tayler-instability driven ones. In these pictures, the magnetic field is typically derived 
from the induction equation, and the field generation is attributed to differential rotation via 
the velocity field of the fluid in a mean-field dynamo [156,189-193]. 


4.3.4. Adaptation in Evolutionary Models 


Despite the various forms and uncertainties of dynamo descriptions, particularly about 
their saturation strengths and resulting viscosities in stably stratified stellar regions, they are 
commonly adopted in stellar evolution models of high-mass stars. To a large extent, they result 
in the favourable outcome that the rotation profile is flattened, e.g., [194], and the core angular 
momentum can be impacted to match the spin rates of compact objects. The observational 
verification of these processes is still challenging. The best way to gain more insights may be 
via asteroseismology of high-mass stars [79,155,159,195—-198]. 


4.4. Fossil Field in the Stably Stratified Envelope 
4.4.1. Magneto-Convection 


Given the vigorous turbulent convection in the cores of main-sequence massive stars, it is 
likely that fossil fields cannot be relaxed in that region. Early analytical works have practically 
neglected this constraint and generally assumed that a large-scale magnetic field could settle 
in a star, permeating through all layers [128]. In fact, even more recent simulations assume a 
self-gravitating non-rotating gas, where the distinct main-sequence structure of a star is not 
explicitly accounted for [114,117]. While this constraint was considered in the simulations 
of [169], and in a phenomenological way in the studies of [124,199], the physical nature of fossil 
fields near convective cores remains largely unexplored. 

Moreover, hot and high-mass stars on the main sequence may have a very thin layer close 
to the surface where some small fraction of the total luminosity is transported via convec- 
tion [200]. This is expected to happen as a result of the iron opacity peak. However, a fossil 
magnetic field may suppress or completely quench convective turbulence in this region. The 
observational verification of this phenomenon might be possible by studying stellar oscillations. 
OB stars in general show evidence for measurable macroturbulent broadening [201]. This 
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also includes known magnetic OB stars, where typical macroturbulent values range from 5 to 
70 km s7! [57,62], which on average is lower than in the total sample of OB stars. The only 
exception is NGC 1624-2 with the strongest measured magnetic field in the sample, where the 
macroturbulence is between 0 and 3 km s ~! [142]. The origin of macroturbulence remains 
debated. Since fossil fields can inhibit convection by the criterion given in Equation (8), it 
seems challenging to explain that, in these cases, macroturbulence could arise from the sub- 
surface region. Instead, it might be more likely that the observed macroturbulence in known 
magnetic OB stars originates from g-mode pulsations from the stellar core [195,202], which 
however could be damped by the magnetic field. The rotational velocities in general are lower 
in magnetic high-mass stars in comparison to those stars where strong surface fields are not 
detected. Since macroturbulence may also correlate with stellar rotation, this might explain the 
lower measured macroturbulent velocities in magnetic high-mass stars. 


4.4.2. Surface Magnetic Braking 


Stars with fossil fields embedded in their radiative envelopes are subject to significant 
angular momentum loss due to magnetic braking (Section 3.4). The spin-down of the stellar 
surface is clearly evidenced from observations [203-205] and MHD models [163,164]. 

The first massive star evolutionary models taking into account this effect were studied 
by [206] using the Geneva stellar evolution code. They found that models with magnetic 
braking strongly depend on the assumptions regarding the internal coupling in the star. To 
further investigate this effect, the authors of [207] also computed Geneva stellar evolution 
models, including magnetic field evolution (in the form of Equation (3), allowing for the surface 
strength of the magnetic field to change as a function of the stellar radius) and magnetic mass- 
loss quenching (see below). The main conclusion is that magnetic braking can rapidly brake 
the surface rotation of the star; however, the magnetic properties and the efficiency of rotational 
braking also change as a function of time. Subsequent modelling approaches using the MESA 
software instrument could investigate a larger parameter space, and hence the dependence on 
initial mass and initial rotation [123]. Some studies also explored surface magnetic braking in 
cases where an internal dynamo-driven magnetic field operates [156,189,191]. 

Another observable outcome of magnetic braking concerns the mixing of chemical ele- 
ments. As the star spins down, the time scale of magnetic braking has to be longer than the time 
scale of chemical mixing to see any surface enrichment of core-produced materials, 14N specif- 
ically. Quantitatively, this depends on the initial magnetic field strength, its time evolution, 
and further assumptions about the angular momentum distribution. Correlations have indeed 
been difficult to identify in observed samples of magnetic stars, e.g., [208-210]. While the direct 
impact of the magnetic field on chemical element transport is not fully clear, stellar evolution 
models incorporating magnetic braking have shown that the spin-down indirectly impacts 
the surface chemical composition. In particular, magnetic models can produce slow-rotating 
nitrogen-enriched stars with high surface gravities, which is not possible with standard evo- 
lutionary models otherwise [124,207]. The existence of such stars is well-documented from 
large-scale spectroscopic studies in the Magellanic Clouds, e.g., [211-214]; however, their 
magnetic characterisation still awaits (see Section 2.2). 


4.4.3. Mass-Loss Quenching 


In pioneering MHD models it was shown that the stellar magnetosphere plays a crucial 
role in the dynamics of the outflow from hot massive stars [44]. Importantly, channelling of 
wind material along magnetic field lines leads to shocks and X-ray emission in the magnetic 
equator [47]. Subsequent analytical and numerical magnetospheric models also evidenced that 
the mass trapped in the magnetosphere eventually falls back onto the stellar surface [163,164, 
215-218], while break-out events may release some plasma [89,219,220]. The net effect is that 
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the total mass outflow is largely reduced by magnetospheric dynamics. The phenomenon is 
also clearly evidenced in terms of the lower derived stellar mass-loss rates [221]. Considering 
the long-term impact of such magnetic mass-loss quenching, stellar evolution models could 
be studied for the first time accounting for reduced mass-loss rates attributed to an evolving 
dipolar magnetic field [222,223]. The evolutionary impact of magnetic mass-loss quenching 
is modest for stars with weaker winds, that is, stars initially less massive than about 20 Mo 
at solar metallicity [207,224]. However, more massive stars can retain close to all of their 
initial mass if a strong magnetic field is present [123,124]. This approach led to subsequent 
applications in the context of heavy-stellar mass black holes, pair-instability supernovae, and 
luminous blue variables [225-227]. 


1 100% 
3 LMC ine 
T 60% 
mA 
F SMC 40% 
Ẹ HD 37022 


oD 
© A 
N 1 1/10 Zo 20% 's 
NX `~ 
<= 
10% 

NGC 1624-2 8% 

6% 

gh 4% 


10 2 4 6 8 10 12 14 
R,/R, (at solar metallicity) 


Figure 3. Equivalency curve between different physical processes leading to a reduction in stellar 
mass-loss rates (colour-coded). The horizontal axis shows the Alfvén radius (in units of the stellar 
radius), quantifying the strength of the magnetic field. It is considered at solar metallicity. Some known 
magnetic massive stars at solar metallicity are included for reference. The vertical axis shows the effect of 
metallicity for a non-magnetic case. For reference, the Large and Small Magellanic Clouds (LMC and SMC, 
respectively) are indicated. The first LIGO detection was inferred at 1/10 of solar metallicity. Adopted 
from Figure 7.2 of [223]. 


The key principle is illustrated in Figure 3. The winds of hot stars are radiatively-driven, 
primarily due to iron-group elements, e.g., [228-231]. Mass-loss rates can be impacted by 
various processes and factors. Considering a reference non-magnetic model at solar metallicity 
(top left corner of the diagram), its mass-loss rate can be decreased by two independent 
processes. Lower metallicity leads to a reduction in M since there are fewer metals to drive 
the wind. A typical metallicity scaling is estimated to be x Z™ with m in the range of 0.70- 
0.85, e.g., [232]. This applies not only to OB star winds but also to Wolf-Rayet type winds 
(which can take away a large amount of mass prior to core collapse) that may be largely 
impacted by metallicity effects [233-235]. Stronger magnetic fields lead to a reduction in M 
since more mass can be trapped in the magnetosphere and channelled back onto the surface. 
This can be quantified with a scaling of ~ 1 — y1 — 1/R, with the closure radius proportional 
to the Alfvén radius Re ~ R, + 0.7(R,4 — R,) for a non-rotating dipolar case [163,164]. The 
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implication, therefore, is that astrophysical phenomena that were thought to be restricted to 
low metallicity environments, due to requiring low mass-loss rates, might be possible at high 
metallicity if the progenitor star is strongly magnetised. 


4.5. Stellar Mergers 


Stellar mergers are exciting astrophysical phenomena. Although these dynamical events 
have not been directly observed yet, there are interesting candidates for post-merger objects, 
e.g., [236-238]. The first works suggesting stellar mergers as the origin of observed magnetism 
in high-mass stars were by [239,240]. Subsequently, the authors of [154,241,242] scrutinised 
whether main-sequence mergers could explain the apparent rejuvenation of stars, developed 
state-of-the-art MHD simulations, and followed the evolution of the merger product. In the 
pioneering MHD model, a seed magnetic field is greatly amplified via the magneto-rotational 
instability during the merger of two stars. This is a very promising channel to explain how 
some stars may become strongly magnetised, and there is support for such a channel for the 
example of compact objects [75,243]. 

The nature of the resulting magnetic field is not immediately clear. The magneto-rotational 
instability creates a primarily toroidal magnetic field, which is likely close to being perpen- 
dicular to the rotation axis, suggesting that there might be some preferential alignment of the 
resulting magnetic geometries. Given the immense computational cost of such simulations, the 
magnetic field evolution could be followed for 5 Alfvén crossing time scales [154]. This is much 
less than the order of few 10-100 Alfvén crossing time scales studied for the relaxation of fossil 
fields [73,114,117,118]. It is thus not yet fully clear if the merger-induced field amplification 
actually does relax into a stable fossil configuration. This might sensitively depend on how fast 
the outer stellar envelope can become stable against convection and differential rotation. 

It was claimed that the models broadly represent the properties of the B-type star t Sco [65, 
115] as a proof of concept [154,242]. As they noted, some observables were not well matched 
with the post-merger models, particularly the observed nitrogen abundance and very slow 
surface rotation of the star. Single-star evolutionary models were studied focusing on the 
reconciliation of these mismatches [244]. They found that order of magnitude changes are 
required in efficiency parameters to speed up magnetic braking and chemical mixing. This 
is indeed unlikely to be compatible with a single-star scenario and could perhaps be better 
explained via a binary channel. Nonetheless, it remains difficult to assess if the special case of 
T Sco could represent other early-type stars in general. For example, the magnetic field of this 
star is much more complex than that of typical OB stars, e.g., [62,66]. 

Observational evidence suggests that main sequence stellar mergers are unlikely to be 
responsible for generating all fossil fields in high-mass stars. Many known magnetic massive 
stars are incompatible with such an evolutionary history; in particular, the existence of some 
young binary systems on stable, nearly circular orbits with one magnetic component, e.g., [245] 
and the doubly magnetic € Lupi system [246] pose a significant challenge. Some statistical 
arguments also point toward a discrepancy to explain the origin of fossil fields solely via merger 
events. The incidence rate of binarity is expected to increase towards higher masses [247]. 
Therefore, the occurrence of stellar merger events is also expected to be mass-dependent, with 
increasing rates towards higher masses [248]. On the other hand, the incidence rate of observed 
magnetism is approximately constant across OBA spectral types and inferred masses (see 
Section 2.2). Moreover, if merger remnants were all strongly magnetised, then certain groups 
of stars (for example, blue stragglers) should show evidence of high magnetic incidence rates. 
This is however not observed [57,249]. 

In summary, the merger of two main-sequence stars represents a possible channel to form 
large-scale fields in high-mass stars. However, it remains disputed whether it could be a general 
explanation for the origin of fossil fields in all high-mass stars. Ample observational evidence 
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(magnetised pre-main sequence stars, randomly distributed obliquity angles, magnetised 
components in close binary systems) points toward the direction that fossil fields must remain 
from earlier stages of the star’s history. Thus, the exact origin of fossil fields in high-mass stars 
is still an unresolved problem. 


5. Overall Picture of Magnetic Field Evolution in a High-Mass Star 


Let us now take a simple, qualitative, and phenomenological overview picture of how 
magnetic fields evolve through stellar evolution, from early to late stages. This concerns various 
spatial and temporal scales, and many details remain highly uncertain. Figure 4 shows a few 
stages of interest. We can assume that the star is a typical high-mass star of around 20 Mo, and 
it evolves through a Red Supergiant phase to core collapse. 


Phase 1: Collapse of the Molecular Cloud 


Molecular clouds are magnetised environments. During their collapse, magnetic fields are 
wound-up in a disk around the forming star. The proto-star accretes material—and perhaps 
magnetic flux—from the surrounding magnetised disk. Dynamo-driven magnetic fields may 
be initiated in situ in developing stellar convective regions and “freeze out” to radiative regions. 
Thus far, the magnetic properties of the forming star are very uncertain. 


Phase 2A: Observably “Non-Magnetic” OB Star 


As a main sequence star forms, its structure is relaxed into a convectively unstable core 
and a stably stratified envelope. The convective core is shown to induce a dynamo action in 
numerical simulations of A- and B-type star models (Section 4.1); however, this phenomenon is 
not yet observationally evidenced. 

The radiative envelope could host a fossil magnetic field without any surface manifestation 
(Section 3.3). If initially differential rotation is present in the system, a dynamo action can also 
be realised. Dynamo-produced magnetic fields are not expected to have any significant surface 
components either. Magnetic fields (either fossil or dynamo) in radiative stellar envelopes 
have not been directly verified via observations. Their inference relies on asteroseismic studies 
or much more indirect and uncertain effects, for example, the overall spin evolution of the 
star. These stars, with a lack of direct detection of surface magnetic fields, are often termed 
“non-magnetic”, although some form of magnetism presumably exists in their interiors. 


Phase 2B: “Magnetic” OB Star 


In a “magnetic” massive star, there is evidence for a strong surface component of the 
magnetic field, as well as co-rotating wind material confined in the magnetosphere (Sections 2.2 
and 4.4.3). The magnetic field is most likely of fossil origin. It must be embedded inside the 
star, satisfying any long-term stability criteria. While the fossil field permeates most or all of 
the stellar envelope, it is unlikely that it could remain stable in the convective core. In fact, 
a core dynamo could well take place in such stars, with a complex interaction at the core 
boundary [169]. The fossil field enforces rigid rotation in those stellar layers where present; 
hence, the source of dynamo action in stably stratified regions is absent. 


Phase 3A: Red Supergiant with a Convective Core 


As the OB star finishes core hydrogen burning, it becomes a Red Supergiant (RSG). In the 
core of the RSG, helium burning will take place. The stability of the core against convection 
depends on the mass, rotation, and metallicity. A convective helium core might induce a 
dynamo action for the example of other convective dynamos. Differential rotation in stably 
stratified regions could also initiate a dynamo cycle. Such internal properties of RSGs remain 
speculative and poorly constrained. 
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Figure 4. Simplified phenomenological picture of magnetic field evolution from early to late stages of 
stellar evolution for a typical high-mass star. The objects are not to scale. The time evolution proceeds 
from the molecular cloud phase towards the compact remnants. The black lines indicate magnetic field 
lines. For stars: Convective burning regions are shown in yellow. Convective non-burning regions are 
shown in red. Radiative regions are shown in blue. The distinct structure of the pre-supernova Red 
Supergiant is greatly simplified. See text for the description of the phases, labelled from 1 to 5. 
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Phase 3B: Red Supergiant with a Convective Core 


RSGs from “magnetic” progenitors could preserve weak (mostly dipolar) fossil fields 
at the stellar surface [207], if those are able to reconnect along the newly developed, large 
convective cells. The convective surface layers are a likely source of dynamo action [250], which 
could manifest in weak and complex magnetic fields as observed, for example, in Betelgeuse 
and some other RSGs, e.g., [27,251]. 


Phase 4: Advanced Shell Burning Stage 


As the star undergoes various burning stages, the structural changes sensitively impact 
the magnetic properties. Fossil fields are essentially unexplored in detail in these stages. The 
convective silicon—oxygen shell is shown to induce a strong dynamo action in numerical 
simulations (Section 4.2). Although the changes in the core are rapid, the stellar surface is 
possibly still characterised by a weak dynamo activity. Models and observations are still scarce 
to draw firm conclusions. 


Phase 5A: Formation of a Typical Neutron Star 


The magnetic field of a neutron star is either generated at around the time of core collapse 
or inherited from the precursor. A generating mechanism is not yet clear, but it is possible that 
neutrino-driven convection or crust formation in the neutron star are key components. 


Phase 5B: Formation of a Magnetar 


Similarly to the uncertainties of how typical neutron stars obtain their magnetic fields, the 
field origin in their subclass of magnetars is also unresolved (Section 2.3). Arguments based 
on magnetic flux conservation suggest that OB stars with fossil fields could be the precursors. 
Dynamo models need to invoke a rapid field amplification to reach an order of 10° G strength 
and saturate at that amplitude. This requires a rapidly-rotating progenitor, which favours 
a weakly magnetised or rotationally weakly coupled core-envelope configuration. In that 
case, the progenitor can preserve a larger core angular momentum reservoir. Alternatively, 
convection or fallback could also power a dynamo mechanism. It is at the forefront of current 
research to solve these open questions. 


6. Conclusions 


Stellar magnetism has been a much studied area of astrophysics. There is now renewed 
interest in this field given the understanding that it might drastically impact all stages of stellar 
evolution and the formation of compact objects. The impact on stellar populations and galactic 
archaeology is still largely unexplored. Major progress has been possible with the advent of 
extensive photometric, spectropolarimetric, and multiwavelength spectroscopic observations. 
It is essential to continue these efforts and survey extragalactic high-mass stars. Direct magnetic 
field detection via the Zeeman effect requires consideration of instrument design and equipment 
of the largest optical telescopes with a high-resolution spectropolarimetric unit. This could help 
study the metallicity dependence of magnetism and answer the puzzling origin of observed 
magnetic fields in high-mass stars. The detected magnetic fields in main sequence high-mass 
stars are understood to be of fossil origin. Their impact on stellar evolution has only just begun 
to be characterised in detailed model calculations. The first large-scale grid of evolutionary 
models, structure models, and isochrones incorporating fossil field effects of massive stars 
has recently been published [124] and made available at: https: / /zenodo.org/record/706976 
6 [252]. Extensive work has focused on dynamo mechanisms in radiative stellar layers. Such 
dynamo models rely on one of three major pillars: the Tayler instability, the magneto-rotational 
instability, and the mean-field theory. Thus far, the verification of these models is difficult; 
however, they may play a crucial role in describing the star’s angular momentum evolution 
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from formation to core collapse. Although dynamos in convective stellar layers are likely to 
exist, they remain almost completely unexplored in several phases of stellar evolution. The 
open questions in this research field motivate further investigations. Magnetic fields can largely 
impact the physics and evolution of massive stars from their formation until core collapse. 
Therefore, it remains a contemporaneous topic in stellar astrophysics. 
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